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ABSTRACT 

Irregular dust grains are subject to radiative torques when irradiated by 
interstellar starlight. It is shown how these radiative torques may be calculated 
using the discrete dipole approximation. Calculations are carried out for one 
irregular grain geometry, and three different grain sizes. It is shown that 
radiative torques can play an important dynamical role in spinup of interstellar 
dust grains, resulting in rotation rates which may exceed even those expected 
from H2 formation on the grain surface. Because the radiative torque on 
an interstellar grain is determined by the overall grain geometry rather than 
merely the condition of the grain surface, the resulting superthermal rotation 
is expected to be quite long-lived. By itself, long-lived superthermal rotation 
would permit grain alignment by normal paramagnetic dissipation on the 
"Davis-Greenstein" timescale tdg- However, radiative torques arising from 
anisotropy of the starlight background can act directly to alter the grain 
alignment on times short compared to tdg- Radiative torques must therefore 
play a central role in the process of interstellar grain alignment. 

The radiative torques depend strongly on the grain size, measured by Oefr, the 
radius of a sphere of equal volume. In diffuse clouds, radiative torques dominate 
the torques due to H2 formation for aes = 0.2/im grains, but are relatively 
unimportant for acs < 0.05/im grains. We argue that this may provide a natural 
explanation for the observation that Oefj ^ 0.1/im grains in diffuse clouds are 
aligned, while there is very little alignment of a^s ^ 0.05/im grains. We show 
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that radiative torques are ineffective at producing superthermal rotation within 
quiescent dark clouds, but can be very effective in star-forming regions such as 
the Ml 7 molecular cloud. 

Subject headings: ISM: Dust, Extinction - Polarization - Scattering 

1. Introduction 

Polarization of starlight by aligned interstellar dust grains was discovered nearly 
half a century ago (Hiltner 1949a,b; Hall 1949; Hall & Mikesell 1949), but the processes 
responsible for the observed alignment remain uncertain. Davis & Greenstein (1951) 
observed that interstellar grains were expected to be rotating rapidly due to Brownian 
motion, and proposed that these spinning grains could be aligned with the local magnetic 
field by paramagnetic dissipation. However, further study of the statistical mechanics of 
grain alignment (Jones & Spitzer 1967; Purcell & Spitzer 1971) raised questions about the 
ability of the interstellar magnetic field to achieve the observed degree of alignment, since 
random gas-grain collisions would tend to oppose the alignment process. 

Purcell (1975, 1979) first recognized that interstellar grains were expected to have 
superthermal rotational velocities as the result of systematic torques, the most important of 
which appeared to be due to the process of H2 formation on the grain surface. As discussed 
by Purcell, superthermal rotation due to torques which are fixed in body coordinates 
can enhance the degree of grain alignment, since the "thermal" torques due to random 
collisions with gas atoms now have little effect on the direction of the grain angular 
momentum, allowing paramagnetic dissipation to inexorably bring the angular momentum 
into alignment with the galactic field. 

The systematic torques considered by Purcell were due to processes taking place at 
the grain surface - H2 formation, photoelectric emission, and inelastic collisions with gas 
atoms - but the grain surface may be altered due to contamination or erosion on relatively 
short time scales. The resulting changes in direction of the systematic torque can cause the 
grain to occasionally undergo short periods when its rotation is "spun down" ; during these 
"crossover" episodes the grain may become disaligncd (Spitzer & McGlynn 1979). Because 
of this disorientation during "crossover" , it was not clear whether ordinary paramagnetic 
dissipation plus superthermal rotation driven by "Purcell torques" can account for the 
observed grain alignment. In addition, dust grains were observed to be aligned in some 
dense molecular regions where "Purcell torques" were expected to be ineffective due to a 
low H/H2 ratio, attenuation of the ultraviolet radiation required for photoelectric emission. 
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and near-equality of gas and grain temperatures. As a result, there has been renewed 
interest in alternatives, including the possibility that grains may be superparamagnetic 
(Jones & Spitzer 1967; Duley 1978; Martin 1995; Goodman & Whittet 1995), or that 
grain alignment is due to gas-grain streaming (Gold 1952; Lazarian 1994, 1995a; Roberge, 
Hanany, & Messinger 1995). Lazarian (1995b) has emphasized the possible importance of 
grain "helicity" , since helical grains can be driven to superthermal rotational velocities, and 
possibly aligned, when exposed to either streaming gas atoms or anisotropic radiation. 

Harwit (1970a,b) suggested that the quantized angular momentum of the photon could 
lead to rapid rotation of a grain following absorption and emission of many photons, and 
proposed that the anisotropy of starlight could result in a tendency for interstellar grains 
to spin with their angular momentum vectors parallel to the Galactic plane. Dolginov 
(1972) observed that interstellar grains might have different absorption and scattering cross 
sections for left- and right-handed circularly polarized light, so that the grain angular 
momentum could be changed if illuminated by unpolarized but anisotropic radiation. 
This effect was further discussed by Dolginov & Mytrophanov (1976) for particles in the 
Rayleigh limit, and by Dolginov & Silant'ev (1976) for larger particles but with refractive 
index close to unity so that the Rayleigh-Gans approximation could be used. Dolginov 
& Mytrophanov noted that this process could lead to both rapid rotation and possible 
alignment. Unfortunately, Dolginov and collaborators were unable to calculate the torques 
on grains with realistic compositions and sizes. 

That irregular interstellar grains should be subject to radiative torques is not surprising. 
We consider two macroscopic examples for illustration. Fig. |l|a shows a four-fold symmetric 
target with square top and bottom, with each of its four rectangular sides divided into 
perfectly absorbing and perfectly reflecting halves. In the geometric optics limit the normal 
component of the radiation pressure force will be twice as large on the reflecting sections 
as on the absorbing sections; as a result, an isotropic radiation field illuminating this target 
will produce a positive torque along axis ai. However, incident radiation which is either 
parallel or antiparallel to ai will not produce any torque on this target. 

Fig. |l|b shows a target obtained by starting with the shape of Fig.|l|a and removing four 
wedges from the top and four from the bottom, with the resulting shape being symmetric 
under reflection through the centroid. Fig.^b is an example of a shape with "helicity". 
From symmetry it is clear that an isotropic radiation field will produce no torque on this 
target: whatever torque is exerted on the "top" half of the grain will be cancelled by an 
opposite torque on the "bottom" half. However, anisotropic illumination can produce a 
torque on this target. If, for example, the surfaces are all perfectly refiecting, then incident 
radiation which is antiparallel to ai will produce a torque parallel to ai, while radiation 
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parallel to ai will produce a torque antiparallel to ai. These two examples show that 
macroscopic objects will be subject to radiative torques unless they are highly symmetric; 
thus irregular targets should generally be subject to radiative torques. Interstellar grains 
are, of course, comparable to or smaller than the wavelength of the illuminating radiation, 
but one does not expect the radiative torques to vanish when target geometries such as Fig. 
ID are reduced to sizes comparable to the wavelength. 

In this paper we discuss the forces and torques on grains of arbitrary shape and 
with sizes which are neither large nor small compared to the wavelength of the incident 
radiation. We show how these forces and torques can be calculated using the discrete dipole 
approximation. Quantitative results are obtained for one particular irregular grain shape. 



The first part of this paper, §3^13) is devoted to development of the theory of scattering 
by irregular particles to enable efficient computation of radiative torques. This is carried out 
within the conceptual (and computational) framework of the discrete dipole approximation. 
In §1^ we report results of extensive computations for one specific irregular grain geometry 
and composition, and three different sizes, agg = 0.2, 0.05, and 0.02/xm (aeg is the radius of 
a sphere of equal volume). 

Readers primarily interested in the implications for grain rotation may elect to skip 
§§@-||, and proceed directly to §1^, where we discuss the effects of starhght and gas drag under 
realistic conditions in interstellar diffuse clouds, and §& where the resulting superthermal 



rotation is evaluated. It is seen in § ^.3] that moderately anisotropic starlight can torque 
Oeff = 0.2/im grains in diffuse clouds up to extremely large rotational velocities, exceeding 
even the superthermal rotation due to H2 formation. Smaller [a^s ~ 0.05/im) grains, on the 
other hand, are only weakly affected by radiative torques. 

We also examine the importance of radiative torques within quiescent dark clouds 
and within active star-forming regions (§ ^.5|) . 



Our results are summarized in §|^. Our principal result is that the torques exerted 
on interstellar grains by background starlight are dynamically very important for grains 
with ttefr ^ 0.1/im. These torques will produce extreme superthermal rotation of interstellar 
grains in both diffuse clouds and star-forming clouds. In a separate paper (Draine & 
Weingartner 1996) we examine the role of these radiative torques in the alignment of 
interstellar grains with the galactic magnetic field. 
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2. The Discrete Dipole Approximation 

We adopt the "discrete dipole approximation" (Draine & Flatau 1994), and represent 
the target by an array of polarizable points, at locations r^, with electric polarizabilities 
aj. The space between the points is vacuum. We consider the response of this array 
to monochromatic illumination, and adopt the usual complex representation for all 
time-dependent variables which are first-order (i.e., E and B fields, and dipole moments 
Pj). Second-order variables (i.e., forces F and torques F) will be taken to be purely real. 
We assume an incident plane wave 

Einc = ^incfl exp(ik-r - iujot) (1) 

Bine = Bi„c,o exp(ik-r - iuot) (2) 
Bjnco = k X Ejjico (3) 

where the unit vector k = k/A;, and k = Uq/c. By suitable choice of complex Ej„c,o we can 
represent general elliptical polarization. 

The dipole at location r^ acquires a dipole moment 

p,(t) = p,(0)e--''* = a,E, (4) 

where E^ is the electric field at r^ due to the incident electromagnetic field and all dipoles 
except dipole j: 

Ej Ej^c J -|- Egcaj (5) 

Eincj = Ej„c,o exp(ik-rj - icuot) (6) 

^ikrji ( / xlzT 7^ r n 1 

E.ea,, = e-^-«* E ^ k'^oi X (P^ X + ^ o ' ' 3r,,(r,rpO - rjpj , (7) 

where r^^ = r^ — r^. Similarly, the magnetic field at r^ may be written (see, e.g., Jackson 
1975). 

-^mCjj ~l~ ^sca,j (8) 

Bj„cj = Bj„c,o exp(ik • r^- - iujot) (9) 
B.caj = e-°* E /c^^(r,, X p,) (l - ^) , (10) 
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3. Electromagnetic Force on a Grain 

The instantaneous force on point dipole j is 

/dp A 



Re (pj-Vj)Re(Ej) + -Re 



dt 



X Re(Bj 



The first term in ( [TT| ) is due to gradients in the local electric field; the second is the Lorentz 
force on the currents associated with the oscillating dipole moments. 



It is convenient to separate the total force Frad = J2j Fj into two terms: 

Frad 



F- + F 

' mc ~ ^ sea ; 



N 

^- = T 

mc / , 

i=i 

N 



Re (p, - V,) Re(E,„,J + ^Re 1^ 



X Re(Bjncj) 



F = V 

■ sea / , 

The time-averaged force is 



Re (pj-Vj) Re(E,ea,i) + -Re 



' dpj 
~dt 



X RefB 



sea,] I 



rad) (f jj2c) ~l~ sea) 

Letting x* denote the complex conjugate of x, it is straightforward to show that 



-Re 



N 



^ ik(p* (0) ■ Ej„c,o) exp (ik ■ ) 



Stt 



(12) 
(13) 

(14) 
(15) 

(16) 
(17) 



where C^xt is the total extinction cross section, evaluated using the optical theorem (Draine 
1988). Thus (Fj„c) is just the average rate at which momentum is removed from the 
incident radiation field. 

Direct evaluation of (F sea) from (|1^), using (|^) and (p^, involves summing over 
A^(A^ — 1) terms, which for ^ 10'^ becomes computationally prohibitive.]^ Alternatively, 
we may evaluate the net rate of momentum transport to infinity by the scattered radiation 
field; let this be denoted l^outj- Overall momentum conservation requires that the 



^ These sums could presumably be evaluated in 0(iV In iV) operations using FFT techniques similar 
to those used for evaluation of the fields Ej (Goodman, Draine, & Flatau 1991), but this has not been 
implemented here. 
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momentum removed from the incident beam either be transferred to the grain or carried off 
by the outgoing scattered radiation: 

— C'ext|Ej„c,Opk = {Frad) + (Fout) ■ (18) 

Stt 

Thus we see that (Font) = —{Fsca)'- Fgca can be thought of as simply the "recoil" of the 
grain due to the scattered radiation. 

The time-averaged rate of transport of momentum by the scattered radiation is 

{Font) = fr^dn^Re{F:^^xBrad) (19) 



- E[P.(0)-n(n-p,(0))]e-^'="- 



dVL n 



(20) 



i=i 

where in ( pUD we have assumed the surface of integration to have radius r ^ A = 27r/fc, 
retaining only the leading terms of 0(r^^) in E and B (see Appendix A). Each scattering 
direction requires evaluation of 0{N) terms; the angular integrations typically require a few 
hundred scattering directions for accurate evaluation. 

We define a dimensionless radiation pressure efficiency vector Qp^ such that the 
radiation pressure force 

(Frad) = QrrrralJ^^^ , (21) 

where the effective target radius agff = (3V"/47r)^/^, where V is the volume of target 
material. The vector Qpj. depends on both the orientation of the grain and the direction of 
propagation and polarization of the incident radiation (e.g., linearly or circularly polarized). 
It is clear that Qpr-k is just the quantity usually described as the "radiation pressure 
efficiency factor" = Qabs + (1 ~ {cos 9))Qsca, where Qats and Qsca are the usual absorption 
and scattering efficiency factors, and (cos(^)) is the mean value of cos^ for the scattered 
radiation, where 6 is the scattering angle (see, e.g., Bohren & Huffman 1983). It is obvious 
that Qpr II k for targets - e.g., spheres - which are symmetric under rotation around an 
axis parallel to k. 



4. Electromagnetic Torque on a Grain 

We choose a coordinate system such that the grain center-of-mass is at r = 0. The 
instantaneous radiative torque on the grain is 

N N 

Trad = E X Fj + E MPj) X Re(E,) . (22) 
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As for the force calculation, we separate the radiative torque into "incident" and "scattered" 
components: 

^rad ^inc ~l~ ^sca j (^"^) 

Tine is the torque exerted on the oscillating dipoles by Ej„c and Bj„c, and Tgca is the torque 
exerted on the oscillating dipoles by the scattered radiation field - i.e., the fields Egca and 
Bsco due to the oscillating dipoles. The time-averaged value of T^nc is straightforward: 



N N 

J2 P* (0) X E,„,,oe^'"^ - X ^ r,- (p* (0) • E,„e,o) e*"^^ 

j=i j=i 



(24) 



As was the case for Fgca, direct evaluation of r^c^ would be time-consuming. By analogy to 
Fscai we appeal to the fact that the time- averaged torque exerted by the scattered field plus 
the time-averaged rate of transport of angular momentum to infinity by this field must sum 
to zero: {Tsca) + (rout) = 0, where Tout is the net rate of transport of angular momentum by 
the scattered radiation field across a fixed surface enclosing the target. In dyadic notation, 
the instantaneous angular momentum flux may be written (Jackson 1975, problem 6.11) 

M = i- (^Re(E)Re(E) + Re(B)Re(B) - x r , (25) 

where I is the unit dyadic. Integrating M over a spherical surface of radius r: 

Tout = / r^^n-[Re(E)Re(E) + Re(B)Re(B)] x n . (26) 

Now, if r ^ A, the field radiated by the oscillating dipoles has |n-Esca| oc r~^, |n-Bsca| oc r~^, 
iTjsca X n| OC r~^, and {'Bgca x n| oc r~^. Thus we need retain only the leading-order terms 
in evaluating each of these quantities. Expanding Egca and B^ca as shown in Appendix A, 
letting r — > oo, 

(r,,„) = -(r„„,) ^-^JdQ Re (5^Vb + S*s\e) , (27) 

Se = Y1 ~ {^■'^j)^ - •Pj(0)exp(-iA;n-rj) , (28) 
j=i L J 

N 

Sb = ii-^PjiO) X rjexp{—ikh-rj) , (29) 

N 

V£;^E[P.(0)-^i(^i•P.(0))]exp(-^M•r,) , (30) 
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N 



\b = y^p,(0) X nexp(— ifen-r,) 
= -n X . 



(31) 
(32) 



Noting that for pure right circularly polarized radiation, the time-averaged angular 
momentum flux is just |Ei„c,opk/87rA; (easily seen by noting that each photon carries 
angular momentum h), it is convenient to define a dimensionless "torque efficiency vector" 
Qr = Qr,mc + Qr,sca such that the time-averaged torque 



IE. 



ottA; ztt 



A 



(33) 



where Urad 
Clearly we have 



Emc,oP/87r is the time-averaged energy density of the incident radiation. 



Q 



r,mc 



Oefr|Einc,o| 



rRe 



N 



N 



^p*(0) X E,„c,oe^''"^ - ^k X ^r,- (p*(0)-E,,c,o) e 

3=1 j=l 



Q 



r,sca 



, (34) 



(35) 



The direction and magnitude of Qr depends on the grain orientation, and on the direction 
of propagation and polarization of the incident radiation. 

It is clear that for right circularly polarized radiation with A ^ Oeff, we will have 

Qr ~ Qabs^- 



5. Target Orientation 

It is necessary to specify the orientation of the target relative to the incident radiation. 
Let the target contain two axes ai and a.2 (with ai-a2 = 0) which are fixed relative to the 

target. Let us further suppose that ai, sl2, and as = ai x sl2 are the principal axes of the 
target's moment of inertia tensor, with moments of inertia /i > /2 > h about these axes. 
The grain orientation is completely specified by the orientations of any two non-parallel 
fixed axes, in particular ai and a.2. 

We define a "scattering frame", defined by unit vectors ei = k, 62 -L ei and 63 = ei x 62. 
Three angles are required to specify the target orientation. The orientation of ai in the 
scattering frame is described by the two angles e [0, tt] and $ G [0, 27i], where 



ai = cos ©ei -I- sin © cos $62 -|- sin © sin $63 



(36) 
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A third angle (3 G [0, 2tt] describes rotations of the target axis sl2 around ai: 

3.2 = — sin 6 cos /9ei + (cos 6 cos /? cos $ — sin /5 sin $)e2 + (cos 9 cos /5 sin $ + sin P cos $)e3 

(37) 

The target orientation in the scattering frame is illustrated in Fig. ^. 

6. Computations for One Shape 

We consider one specific grain shape. This shape, shown in Fig. ^ is not intended to 
be reahstic, but is chosen to be (1) easily described, (2) asymmetric, and (3) well-suited 
to representation by a cubic array of point dipoles. The shape chosen is an assembly of 
thirteen identical cubes. The coordinates of these cubes (in units of cube width) are listed 
in Table |l|. The first eight cubes are arranged to form a single larger cube; the remaining 
five cubes are attached to the central eight cubes. 

Let the eigenvalues h > I2 > I3 of the moment of inertia tensor be written 

Ij = '^ajMals = ^ajpals ; (38) 

where p is the solid density, and aj are geometric factors. A sphere has 
aj = 1; a 2:2:1 rectangular solid ("brick") has ai = {5/3){n/3f/^ = 1.719, 
«2 = tts = (25/24)(7r/3)^/^ = 1.074. The principal axes a^- and factors aj for our 
adopted grain geometry are given in Table ^ 

Let Tgr be the grain temperature. Purcell (1979) has shown that if the grain is spinning 
with angular momentum J, and ^ IikTgr, then viscoelastic damping and Barnett effect 
damping will rapidly bring the grain's principal axis ai into alignment with J. In our 
discussion below we will therefore assume that ai is at all times perfectly aligned with the 
angular momentum J of the grain. The grain rotation period will be short compared to 
all other dynamical time scales, so it is appropriate to average the forces and torques over 
the grain rotation angle /3, and henceforth we will discuss only the "/3-averaged" efficiency 
vectors Qpr(6,<l>) and Qr(6, $). 

We have carried out scattering calculations for grains with the shape shown in Fig. 
^, the refractive index of "astronomical silicate" (Draine & Lee 1984)| and effective radii 
Oeff = 0.2, 0.05, and 0.02/im. The computations were carried out using the discrete dipole 



^ This dielectric function may be obtained by anonymous ftp to astro.princeton.edu, file 
'draine /dust / eps . Sil" . 
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approximation program DDSCAT.Sa^. The targets were represented by arrays of = 6656, 
22464, and 53248 dipoles; in each case, was taken to be large enough to satisfy the vahdity 
criterion \m\kd < 1 (Draine & Flatau 1994), where m{uj) is the complex refractive index, 
and d = {in/SNy^^aes is the interdipole spacing. Lattice dispersion relation theory (Draine 
& Goodman 1993) was used to determine the optimal choice of dipole polarizabilities 
at each wavelength. The computations employed FFT techniques (Goodman, Draine & 
Flatau 1991) and a stabilized version of the bi-conjugate-gradients method (Flatau 1996); 
nevertheless the computational demands were substantial, because of the need to compute 
scattering for many different wavelengths and orientations. For example, determining the 
scattering properties for one orientation and two polarization states for the Oefr = 0.2/im 
target at a single wavelength A = 0.4/xm, with the target represented by = 22464 dipoles, 
required 720 cpu-sec on a Sun-10 computer with a 125MHz Hypersparc cpu. For the 
aeff = 0.2yum grain, computations were carried out for 26 distinct wavelengths (from O.l/xm 
to 20/im) and 285 orientations (15 values of (3 for each of 19 values of 0), requiring ~ 2000 
cpu-hr. 

Consider a grain spinning around axis ai, with ai ± k. In Fig. ^we show 

n ^ a..(E±ai) + a.,(E || aQ 
27ra^ff 

_ CUE ± aQ - CU^ II aQ 

Vpoi = 2 l^Uj 

Trais 

where Cext is the extinction cross section averaged over grain rotation angle (3. We see that 
the polarization cross section for the aes = 0.2/im grain peaks near A = 0.5/im. Grains of 
approximately this size would therefore be suitable for producing the bulk of the observed 
interstellar polarization, which tends to peak near A ~ 0.55/im. The = 0.05/im grain, 
on the other hand, is most effective at polarizing near A = 0.2/im (see Fig. §); the observed 
wavelength-dependence of interstellar polarization indicates that grains of this size are not 
appreciably aligned in the interstellar medium (Kim & Martin 1994b). 

We can restrict our calculations to $ = 0; for other values of the efficiency vectors 
Q may be obtained by a simple rotation: 

Q(0, $) = Q(0, 0)-eiei + Q(0, 0)-e2(e2 cos $ + eg sin $) + Q(0, 0)-e3(e3 cos $ - 62 sin $) . 

(41) 



The Fortran program DDSCAT.5a is available from B.T. Draine and P.J. Flatau; contact 
draine@astro.princeton.edu for further information. 
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Let {Qpr, Q,pr) and (Qr, Qr) be the radiation pressure and radiation torque efficiency 
vectors for polarization states 1 and 2 (these can be be any two orthogonal polarization 
states, e.g., two linear polarization states, or left- and right-handed circular polarization 
states). For unpolarized incident radiation (i.e., equal amounts of two orthogonal 
polarization states) we have efficiency vectors 

= liQlr + Qlr) (42) 

Qr = ^(Qr + Qr) • (43) 

Henceforth in this paper we will be concerned only with unpolarized radiation, and we will 
omit the superscript. In Fig. ^ we show our results for the component of Qpr(0j 0) along 
k, as a function of A, for several values of G. This is, of course, just the usual "radiation 
pressure efficiency factor" [Qabs + (1 ~ {cos Q))Q sea]- We see the expected tendency for 
Qpr-k to be small when A ^ ttes, and to be of order unity when 27raefr/A ^ 1 (i.e., A ^ 1/xm 
for Qes = 0.2/im, and A ^ 0.3/im for a^s = 0.05/im. 

In Fig. we show the three components of the radiation torque efficiency vector Qp 
for the Oeff = 0.2/im target oriented at = 30°. For A > 0.1585/im results are shown for 
three different dipole arrays, with N = 6656, 22464, and 53248 dipoles, with interdipole 
separations d = 171, 114, and 85.7 A. At A = 0.1585/im we have m = 2.153 -|- 0.263i, so that 
for the three dipole arrays at this wavelength we have \m\kd = 1.47, 0.98, and 0.74. We 
see that the DDA calculations appear to have converged when \m\kd < 1. At the shortest 
wavelength A = O.l/im we have m = 1.545 + 0.956i, and for the = 53248 array we have 
\m\kd = 0.98, marginally satisfying our validity criterion. 

Figs. H and |^ show, as a function of A, the component of Qr along the principal axis 
ai about which we have assumed the grains to be spinning, for aes = 0.2yum and 0.05/im 
grains. This component of Qr will tend to spin up the grain if Qr-ai/J-ai > 0, or to spin 
it down if Qr-ai/J-ai < 0. 



7. Interstellar Grains 

7.1. Radiation Fields 

The average interstellar radiation field ("ISRF") spectrum in the solar neighborhood 
has been estimated to be (Mezger, Mathis, & Panagia 1982; Mathis, Mezger, & Panagia 
1983) 

Amisrf,a = U{\) for A < 2460 A (44) 
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E ^iBx{T^) for A > 2460 A (45) 



c i=l 



where the ultraviolet function 



f/(A) = for A < 912 A or A > 2460 A 

= 1.287 X 10"^(A/;um)^-^^^^ ergs cm^^ for 912 - 1100 A 

= 6.825 X 10~^^(A//im) ergs cm'^ for 1100 - 1340 A 

= 2.373 X 10"^^(A//im)-°-^^^^ ergs cm"^ for 1340 - 2460 A . (46) 

The three blackbody components are indicated in Table ^, where we also provide values of 
u and A for each of the four spectral components, where 

u = J UxdX , (47) 

A^^^ . (48) 

The total starlight energy density is uisrf = 8.64 x 10^^^ ergscm"^. We expect appreciable 
anisotropy in the typical interstellar radiation field: generally, there will be more starlight 
from the direction of the Galactic center, the nearest bright star or association may be 
important, and nearby dust may attenuate the starlight over one part of the sky. We 
define an anisotropy parameter 7, such that the net Poynting flux is •yUrad.c, where Urad is 
the total energy density of radiation. We assume 7 ~ 0.1 to be a characteristic value. In 
our discussions below we will represent the radiation field as an isotropic component with 
energy density (1 — 'y)urad plus a unidirectional component with energy density •yurad- 



7.2. Radiative Forces and Torques due to Starlight 

7.2.1. Unidirectional Radiation 

For a grain subject to a unidirectional radiation field with energy density 'jUrad, the 
force and torque on the grain are 

Frad = 7ra^ff7Mrad(Qpr(6, $)) (49) 

Frad = nals lUrad ^ {QriQ,^)) (50) 
where () now denotes spectral averaging: 
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In Table ^ we present values of ai- (Qp^) and ai- (Qpr) for (6, $) = (0, 0), with the 
spectral averages evaluated for each of the ux in Table 0, for three grain sizes, Oefr = 0.02, 
0.05, and 0.2yum. Note that the torque efficiencies (Qr) become very small as the grain size 
is decreased below ~ O.l/im. 

In Figs. |lO|-[l^ we show the three vector components of (Qr(6, 0)), for the = 0.2 
and aeff = 0.05/im grains; results are again shown for each of the ux in Table ^ 



7.2.2. Isotropic Radiation 

Consider a grain spinning around axis ai. An isotropic radiation field will result in a 
force Frad and torque Trad on the grain which are each fixed in body coordinates. Because 
of the grain rotation, in inertial coordinates the only component of the force or torque 
which will not average to zero will be the component parallel to the rotation axis ai: 

Frad = 7ralsil-l)Urad{Q;r)^l (52) 
Trad = 7Talsil-l)Urad^{Q'n^l (53) 

where 

(Qi-) = - fsinerfe ai-(Q(e,0)) (54) 
2 Jo 

is the effective efficiency vector for a grain spinning around ai and exposed to isotropic 
radiation. Values of {Qp°) and {Qf°) are given in Table |. 



7.3. Rotational Damping 

For a grain in neutral gas of hydrogen density nn = n{}i) + 2ri(H2) and temperature 
T, the gas drag torque on a grain with angular velocity u around axis ai may be written 

2 

Fdrag,gas = - -^^h ( 1.2) (STTmn/cT) ^/^OggC^ai . (55) 

where mn is the mass of an H atom. If all impinging atoms "stick" and then leave with 
negligible velocity relative to the local (moving) surface, then 6 = 1 for a sphere, and 
5 = 2(77/3)^/^ = 2.01 for a 2:2:1 brick; in general we expect 5 ~ ai. The factor 1.2 in 
eq. ( |55D allows for the effects of helium with nne = 0.1nH.Q 



We have assumed the H to be fuUy atomic; in fully molecular gas the factor 1.2 in eq. (|5^) should be 
replaced by 2-^/2 ^ o.2 = .907. 
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The rotational damping time is 

where ps = p/3gcm~^, a_5 = acfr/10~^cm, and T2 = T/IO^K. 

In addition to coUisions with gas atoms, there is rotational damping associated with 
absorption and emission of photons by the grain (Purcell & Spitzer 1971; Roberge, deGraff 
& Flaherty 1993). If we assume that the grain is heated by starlight to a temperature T^, 
and that Qabs oc A"'^ at wavelengths A ^ O.lhc/kTd, then the damping time due to thermal 
emission may be written 

8ai(/j + 3)C(/9 + 4) palsjkT.y 
= 1.60 X lOV ai P3 al, (^)' (^) (58) 

V18K/ V Urad / {Qabs) 

where ({x) denotes the Riemann (^-function, 

{Qabs) = / UxQabsWdX , (59) 

Urad 

and we have assumed /? = 2 in eq. (|58D , as expected for simple models (Draine & Lee 1984) 
and as appears to be required by the far-infrared emission from dust in diffuse clouds (Draine 
1994). Q In eq. ( |58D it is assumed that the radiated photons have angular momentum h 
relative to the grain center of mass; this will be true for radii a <^ hc/kTd = 800(18K/Trf)pm. 
There is additional damping associated with absorption of starlight photons, but it is 
smaller than the damping due to thermal emission by a factor ~ T^/Trad ~ 1/500 where 
Trad is the color temperature of the radiation responsible for heating the grain (Purcell 
1979), and hence may be neglected. Values of (Qabs) are given in Table ^ The rotational 
damping time Tdrag is given by 

'^drag ^drag,gas ~^ ^drag,em " (60) 



8. Superthermal Rotation 



^ The damping time Tdrag,em is insensitive to the precise vahie of /?: the factor (/? + 3)C(/5 + 4)/C(/3 + 3) 
decreases by only 22% if (3 is reduced from (3 — 2 to \. 
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8.1. H2 Formation 

The thermal rotation rate for a grain may be written 

2 15 kT 

oJt = o 5- , (61) 

assuming rotation around ai with kinetic energy kT/2. Purcell (1979) pointed out that H2 
formation on the grain surface, taking place only at randomly- located "active sites" , would 
result in a large torque. The component of the torque perpendicular to the rotation axis 
ai would average to zero; for a 2:2:1 rectangular solid Purcell estimated the component 
parallel to the rotation axis 

TH.-ai = I (lY' f ni}i)i2EkTy/^als I pit) (62) 



3 V3 

where ?t,(H) is the density of H atoms in the gas, is the surface area per H2 formation 
site on the grain surface, / is the fraction of arriving H atoms which depart as H2, -E 
is the kinetic energy of the departing H2 molecules, and p{t) is a random variable with 
time averages {p{t)) = 0, {p{t)p(t + r)) = e""^/*", where to is the "lifetime" of a surface 
recombination site. 

If the only torques acting on the grain are those due to gas drag and H2 formation, the 
grain will attain a rotation rate uJn^; the mean kinetic energy will exceed 0.5/cT by a factor 
(Purcell 1979) 



^HaN^ _ 5ai/^ fn{B.)Y E pa^^l^ ( r^rag y ( to 



Ut J 216(^2 y riH / kT mn \rdrag,gas J \to+rdrag, 



(63) 



8.2. Radiative Torques 



If the grain is subject to a steady radiative torque 

A 



rad — '^(^cS^rad 



27r 



l-7)(Qn+7(Qr(e))-ai 



ai 



(64) 



along the rotation axis ai then, ignoring other sources of rotational excitation by the gas, 
we equate eqs. (p3| ) and ([0^ ) to find that the grain will rotate around ai with an angular 
velocity 



5A 



kT 



1/2 



^rad 



f Urad \ 

SSalfi ySirmii J XrinkT J 



;i-7)(Qn+7(Qr)-ai 



Td, 



rag 



Jd 



(65) 



rag, gas ^ 
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We see that, in addition to (Qr°) (Qr)'^i, uj^ad is determined by the ratio of radiation 
pressure to gas pressure (urad/nukT) and the anisotropy factor 7. The rotational kinetic 
energy will exceed 0.5kT by a factor 



\ UJt J 



192(52 KnnkT) \ nin 



Tdr 



4.72xl09^P3a~5 



/ Urad y 



;i-7)(Qr)+7(Qr)-ai 

"(l-7)(Qn + 7(Qr)-a; 




rag, gas ^ 
2 

1 



(66) 



^T~drag,gas J 
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8.3. Diffuse Clouds 

In a diffuse cloud where the hydrogen is predominantly atomic, H2 recombination on 
the grain surface will result in very rapid grain rotation. Values of [ujh^/ujt)'^ are given in 
Table ^ for three grain sizes, and three possible values of the surface recombination site 
lifetime to. 

A typical diffuse cloud may have nuT = 3000 cm~^K; the average interstellar radiation 
field (ISRF) then corresponds to Urad/nukT = 2.09. We adopt an anisotropy factor 7 = 0.1 
as representative. 

As seen from Figs. |1^ and |ll], the spectrum-averaged (Qr(0))-ai varies considerably 
with 0, the angle between ai and the radiation flux. Values of {urad/^rY are given in 
Table |^ for = 0° and 60°, and for three grain sizes. 

For the agfj = 0.2yum grain at = 0°, we find that radiative torques result in extreme 
superthermal rotation with {urad/^rY ~ 9 x 10"^, greater than the rotation due to H2 
formation on the grain surface if / ^ 35 A. The acs = 0.05/im grain, on the other hand, has 
a radiatively-driven rotation rate only slightly in excess of thermal, with {uOrad/^TY' = 18 
(for = 0°), far smaller than the rotation velocity expected from H2 formation. The 
ttes = 0.02/im grain is essentially unaffected by radiative torques: the rotational kinetic 
energy resulting from these torques is small compared to kT. We therefore see that - 
at least for the particular shape considered here - radiative torques will drive the larger 
(oeflf ^ 0.1/im) grains to highly superthermal rotation rates, while having only a slight effect 
on the smaller (oeff ^ 0.05/im) grains. 

It has been known for some time that the observed wavelength dependence of 
polarization requires that if the small grains and large grains have similar shapes, then the 
degree of alignment declines rapidly for ^ O.lpm (Mathis 1979, 1986; Kim & Martin 
1994a,b). This was somewhat unexpected: the standard Davis-Greenstein alignment time 
varies as tdg oc a^, so one might naively expect that small grains would exhibit a higher 
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degree of alignment than large grains. 

Several ideas have been put forward to explain the low degree of alignment of small 
grains. Mathis (1986) proposed that a grain would not be aligned unless it contained at 
least one superparamagnetic inclusion, and that small grains were statistically unlikely to 
contain such an inclusion. Mathis obtained a good fit to the average polarization curve by 
assuming that a grain with a^s = 0.08/im had a 50% probability of containing one or more 
inclusions. Lazarian (1994) argued that if grain alignment was due to the Gold mechanism 
(gas-grain streaming) driven by Alfven waves, then only the larger grains {aes ~ O.l/zm) 
were sufficiently inertial to be aligned. Lazarian (1995c) argued that H2 formation on fractal 
grains would only be able to drive grains with a^s ~ 0.02yum to superthermal rotation. 

As pointed out by Purcell (1975, 1979), superthermal rotation contributes to grain 
alignment by rendering unimportant the disaligning effects of random collisions with 
gas atoms. Purcell considered three processes which could drive grains to superthermal 
rotation: H2 formation on the grain surface, photoelectric emission, and variations in the 
effective accomodation coefficient over the grain surface; of the three, H2 formation appears 
to be the most powerful. As discussed in §8.1| , the effects of this torque depend on the 
"correlation time" to for the torques driving superthermal rotation. If the only torques are 
due to H2 formation and collisions with gas atoms, the grain will undergo spindown and 
"crossover events" with a frequency t~^, where the time ~ T^itoTdragY^'^ (Purcell 1979). 
If tz -C Tug, then grain alignment by paramagnetic dissipation will be suppressed, whereas 
if tz ^ tdg; then paramagnetic dissipation will result in a high degree of alignment. The 
three processes discussed by Purcell depend upon surface properties; if the grain surface 
is altered on a relatively short time scale, then substantial grain alignment will not occur 
unless the grains are superparamagnetic (in which case t^g is short). 

Here we propose another mechanism for selective alignment of the larger grains. The 
radiative torques on the grain depend upon the global geometry of the grain, and will 
therefore have a correlation time of order the grain lifetime, ~ lO^yr. For the Oeff = 0.2/^m 
grains considered here, ir^adP ^ |rH2p; and we may expect long-lived superthermal 
rotation. By itself, long-lived superthermal rotation would lead to grain alignment by 
paramagnetic dissipation on the Davis-Greenstein timescale tdg, as expected for "Purcell 
torques" which are fixed in body coordinates. However, when an anisotropic radiation field 
is present, the radiative torque components perpendicular to ai may act directly to change 
the grain alignment on time scales short compared to tdq; it appears that this may be the 
mechanism responsible for the observed alignment of aes ~ O.l/im dust grains with the 
interstellar magnetic field (Draine & Weingartner 1996). 

For the a^s = 0.05/im grains, on the other hand, the radiative torques are weak 
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compared to the Purcell torques due to H2 formation. If the surface site hfetime to ^ tdg, 
there will be frequent crossover events, with a resulting low degree of grain alignment. 
This may explain the observed preferential alignment of the larger aes ~ O.lyum interstellar 
grains, and the minimal alignment of ^ 0.05/im grains. 

The dynamics of grain alignment including radiative torques will be examined in more 
detail in a later paper (Draine & Weingartner 1996). 

8.4. Dark Clouds 

The alignment of grains in dense clouds has been problematic. Recent studies of 
polarization of background stars indicate a low degree of grain alignment in the central 
regions of some dark clouds (Goodman et al. 1995). This has been attributed to lack of 
superthermal rotation in regions where the hydrogen is predominantly molecular, where the 
grains are shielded from ultraviolet photons capable of producing photoelectric emission, 
and where the gas and grain temperatures tend to be similar. However, observations 
of polarized FIR emission from M17, M42 (Orion), and W3 (Hildebrand et al. 1995; 
Hildebrand 1996) shows that the grains responsible for the FIR emission in these dark 
clouds are significantly aligned, and it has not been clear why grains in some dense 
molecular regions are aligned, while in others they are not. 

The dark cloud L1755 observed by Goodman et al. (1995) has rin ^ lO'^cm"'^ and 
T !^ 20K, with a peak extinction through the cloud Ay ^ 8. We suppose the cloud to be 
externally illuminated by the average interstellar radiation field. At the cloud surface we 
have Urad ~ 0.5'Uisrf and 7 = 0.5, and at a depth A ^ 2mag from the cloud surface, we 
estimate^ Urad ~ O.OTmisrf, 7 ~ 0.7; the radiation will of course be substantially reddened, 
but we neglect this here. Table ^ lists values of {uradl^xY for grains at this depth. We see 
that the weakened radiation field, and the higher gas density, make the radiative torques 
relatively unimportant in this environment. Near the cloud surface, of course, conditions 
are closer to the diffuse cloud conditions of Table ^, and radiative torques will maintain 
superthermal rotation. Thus we expect grain alignment near the cloud surface, but not at 
depths A ^ 2mag. 



^ If the dust has albedo ~ 0.5 and is highly forward-scattering, then Urad ~ 0.5urarf£^2(0.5v4/1.086), 
7 ~ i?3(0.5A/1.086)/i?2(0. 5^/1. 086), where E2 and E3 are exponential integrals. For ^4 = 2, we would have 
« .0742UISRF, 7 « 0.739. 
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8.5. Star-Forming Clouds 



In star-forming clouds, there are obviously sources of starlight within the clouds. 
To estimate the characteristic starlight energy density to which the FIR-emitting dust is 
exposed, we note that for grains with far-infrared emissivity Qabs oc A~^, the power radiated 
by the grain P oc T|, where is the grain temperature. If we assume that the grain is 
heated by starlight with an energy density Urad, then, since the average interstellar radiation 
field MisRF heats grains to ^ 18K, we find 



As an example, we consider one point on the M17 molecular cloud: {a, 6) = 
(18'^17'"40*, — 15°16'). The lOO/xm emission at this point has a linear polarization of 
~ 5% (Hildebrand et al. 1995), indicating that the grains are appreciably aligned. The 
emission from the dust has a 50//m/100/im color temperature of ~ 85K (Gatley et al. 1979); 
for a emissivity this color temperature corresponds to a dust temperature Td ^ 45K. 
If heated radiatively, these grains must therefore be exposed to a radiation fieldQ with 
Urad/uisRF ~ (45/18)^ = 240 if the radiation heating the grains has a spectrum similar 
to interstellar starlight ("ISRF"). While the spectrum no doubt differs from the ISRF 
spectrum, it is not obvious how: the embedded massive stars will tend to be bluer, but 
the spectrum will be reddened by propagation through the dust. For the present purposes, 
we will assume the ISRF spectrum, with an anisotropy factor 7 = 0.3. The gas in M17 
has riu ~ 10^ cm~^ (cf. Stutzki & Glisten 1990), and presumably has T ^ 45K, close to 
the grain temperature. We use eq.(|^) to obtain the Urad values given in Table ^ We 
see that a^s = 0.2yum grains will have {urad/^TY ~ 2 x 10^ (for = 0°). Evidently the 
radiation field within a star-forming region like M17 is sufficiently intense to drive extreme 
superthermal rotation of the a^fi ^ O.l/im grains. In addition to producing superthermal 
rotation, the components of the radiative torques perpendicular to ai may bring about the 
observed grain alignment (Draine & Weingartner 1996). 



The principal results in this paper are as follows: 



^ For Td = 45K, the 50/um optical depth at this position is r(50/ini) w 0.08, so that the grains must be 
heated by radiation at A <C 50/im. 




(68) 



9. Summary 
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1. We show how forces and torques on an irregular grain may be calculated using the 
discrete dipole approximation. 

2. We report numerical results (in Fig. |^ and Table ^ for k-Qp^,(A), the component of 
the radiation pressure efficiency vector parallel to the radiative flux for one particular 
irregular grain geometry, where the grain is assumed to be spinning around its 
principal axis ai of largest moment of inertia. 

3. We report numerical results for the "radiation torque efficiency vector" Qr(A), (in 
Figs. H - ^) for one irregular grain geometry. The torque depends upon the angle 
between the incident flux and the grain rotation axis ai. 

4. In interstellar diffuse clouds, radiative torques will drive a^s ^ 0.1/im grains to 
extreme superthermal rotation. 

5. Radiative torques will dominate the torques due to H2 formation for a^s ~ O.lyum 
interstellar grains. For these grains, the superthermal rotation is expected to be 
very long-lived, since it depends upon global properties of the grain, rather than 
the relatively short-lived surface properties which determine superthermal rotation 
due to H2 formation, photoelectric emission, or variations in surface accomodation 
coefficient. 

6. In an isotropic radiation field, the long-lived superthermal rotation of aes ~ O.l/xm 
grains due to radiative torques would facilitate alignment of these grains with 
the Galactic magnetic field by the Davis- Greenstein mechanism of paramagnetic 
relaxation. 

7. Since the rotation of smaller (a ^ 0.05/im) grains is not dominated by radiative 
torques, alignment of these grains with the Galactic magnetic field would be limited 
due to random variations in the torque associated with H2 formation on the grain 
surface. This may explain why aes ~ O.l/xm grains in diffuse clouds are aligned, while 
there is evidently minimal alignment of the a^s ~ 0.05/im grains which dominate the 
extreme ultraviolet extinction. 

8. Radiative torques appear to be able to drive grains to superthermal rotation in 
star-forming regions, such as M17, where the ratio of anisotropic radiation to gas 
pressure is relatively high, thereby potentially explaining the observed alignment of 
far- infrared emitting dust in Ml 7 and other star- forming regions with warm dust. In 
cold dark clouds, on the other hand, radiative torques are unable to drive the grains 
to superthermal rotation, consistent with the observed lack of aligned grains deep in 
quiescent dark clouds. 
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9. In addition to the torque component parallel to the grain rotation axis ai (which 
drives superthermal rotation) , anisotropy of the radiation field can result in radiative 
torque components perpendicular to the grain angular momentum; these components 
can afi^ect grain alignment directly. 
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A. Appendix: Radiation Field due to Oscillating Dipoles 



Consider an array of oscillating point electric dipoles pj, at locations r^. At large 
distances r the fields due to these dipoles may be written 
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Bsca X n = — e*'^" J2 [Pj - n(n-Pi)] exp{-tkh-rj) + 0{r-^) (A6) 
^ i=i 

These expressions, substituted into eg. (PBD , lead to the result ( [27| ) for the time-averaged 
torque (Tout)- 
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Table 1. Target Geometry: Coordinates of Constituent Blocks 
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Table 2. Target Principal Axes and factors aj 



j {^j)x {^j)y {^j)z 



1 0.4488 0.4357 0.7802 1.745 

2 0.6673 -0.7441 0.0317 1.610 

3 0.5944 0.5064 -0.6247 0.876 
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Table 3. Interstellar Radiation Field Components 
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Table 4. Silicate Grain Properties 
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Wpr /ISRF 


1.1 X lO""' 


—5 25 X 10""^ 


—7 86 X 10^3 

1 • L/ /\ X \J 




2 64 X 1 0~^ 


2 52 X 1 0"^ 


5 34 X 10"^ 

0.(J^ /n XL/ 


//Oiso\ 
Wr /7500K 


4.43 X 10~^ 


2 04 X 10-6 


5 31 X 1 0-^ 

(J.(JX /\ XL/ 


//Oiso\ 
Wr /4000K 


1 41 X 10~^ 


1 89 X 10-^ 

A. .\jU /n XL/ 


1 06 X 10"^ 

X . L/L/ /N XL/ 


//niso\ 
Wr /3000K 


Q 03 X 10~^ 


6 48 X 10-^ 

W.^tJ /N XL/ 


2 99 X 10-5 

, tj tj /\. XL/ 


//Oiso\ 
Wr /ISRF 


2 84 X 10~^ 


2 75 X 10-5 

^ . \ \J /S XL/ 


1 99 X 10-^ 

J., kJkJ y\ XL/ 


K- \y4pr/uv {'^ — ^ ) 


1 141 


1 sfin 

1 .ouu 


1 .400 


\^pr/7500K I.'-' — ^ j 


u.uzuo 


n 1 84 

U. io4 


1117 


ir Id \ (C^ — (\°\ 


u.uuyoo 


U.U404 


U.OOi 


ij. / r\ \ ICS 

*^-\HprimooK ya — \J ) 


n nn70/i 
U.UU ( Z4 


n n9c;7 
u.uzo / 


n /I c;7 


k-(Qpr)lSRF (6 = 0°) 


0.105 


0.226 


0.784 


ai-(Qr)«. (e = 0°) 


2.54 X 10-=^ 


2.94 X 10-2 


7.01 X 10-2 


ai-(Qr>«. (e = 60°) 


3.55 X 10-3 


1.67 X 10-2 


9.06 X 10-2 


ar(Qr)7oooK (6 = 0°) 


6.01 X 10-6 


2.00 X 10-^ 


7.45 X 10-2 


ar(Qr)7000K (6 = 60°) 


1.22 X 10-5 


2.16 X 10-^ 


4.38 X 10-3 


ai-(Qr)4oooK (0 = 0°) 


1.60 X 10-*^ 


1.88 X 10-5 


7.78 X 10-3 


ai-(Qr)4oooK (0 = 60°) 


3.88 X 10-6 


2.90 X 10-5 


4.39 X 10-3 


ai-(Qr)3oooK (0 = 0°) 


9.01 X 10-^ 


7.36 X 10-6 


2.07 X 10-3 


ai-(Qr)3000K (0 = 60°) 


2.27 X 10-6 


1.52 X 10-5 


2.07 X 10-3 


ai-(Qr)isRF (0 = 0°) 


2.93 X 10-5 


3.59 X 10-^ 


1.69 X 10-2 


ar(Qr)iSRF (0 = 60°) 


4.26 X 10-5 


2.33 X 10-^ 


4.35 X 10-3 
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Table 5. Silicate Grains in Diffuse Clouds 





aeflf = 0.02/im 


ttes = 0.05/im 


ttefr = 0.2yum 


Tdrag,gas (j/^^ 


1.53 X 10^ 


3.81 X 10^ 


1.53 X 10^ 


'^drag,em (j/^^ 


3.12 X 10^ 


2.68 X 10^ 


8.18 X 10*^ 


Tdrag (yr) 


1.03 X 10^ 


3.34 X 10^ 


1.50 X 10^ 


ut (rad s~"^) 


7.01 X 10^ 


7.10 X 10^ 


2.22 X 10^ 


{i^m/uJry for to = lO^yr 


2.10 X 10^ 


4.16 X 10^ 


5.65 X 10^ 


{uJm/uJTf for to = lOVr 


3.86 X 10^ 


1.35 X 10^ 


3.62 X 10^ 


{uJm/uJTf for to = lOVr 


4.22 X 10^ 


1.75 X 10^ 


7.86 X 10^ 


(cOrad/iOTy for 9 = 0° 


3.56 X 10-2 


18.4 


8.73 X 10^ 


{uJrad/'^Ty for e = 60° 


5.49 X 10~2 


11.5 


1.00 X 10^ 



Note. — We assume the grain properties of Table H, p = 3 gem ^, ai = 1.745, 6 = 2, 
nn = 30cm-3, T = lOOK, Ta = 18K, / = 1/3, n(H)/nH = 1, Eu^ = 0.2eV, / = 10 A, 
Urad = MisRF, and 7 = 0.1. 
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Table 6. Silicate Grains in Dark Clouds 





ttcfj = 0.02/im 


ttcff = 0.05/im 


ttcg = 0.2yum 


Tdrag, gas (jf^^ 


1.02 X 10^ 


2.56 X 10^ 


1.02 X 10^ 


'Tdrag,em iy^^ 


9.63 X 10^ 


8.29 X 10^ 


2.52 X 10^ 


Tdrag (jr) 


1.02 X 10^ 


2.56 X 10^ 


1.02 X 10^ 


ojt (rad s""*^) 


3.14 X 10^ 


3.17 X 10^ 


9.92 X 10^ 


{uJm/i^Tf for to = 10 V 


0.466 


1.15 


4.27 


{uJn2/^Tf for to ~ lOVr 


0.470 


1.17 


4.70 


{uJrad/i^Tf for e = 0° 


1.30 X 10-*^ 


9.63 X 10"^ 


4.04 


{Urad/^T? for e = 60° 


5.39 X 10-*^ 


2.10 X 10-^ 


0.275 



Note. — We assume the grain parameters of Table ^ p = 3 gem ^, ai = 1.745, 6 = 2, 
nn = lO^cm-3, T = 20K, / = 1/3, n(H)/nH = .01, E^^ = 0.2eV, / = 10 A, Td = lOK, 
Urad = •07misrf, and 7 = 0.7. 
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Table 7. Silicate Grains in Star-Forming Clouds 





tteff = 0.02/im 


Oefi = 0.05/im 


ttcff = 0.2yum 


^drag,gas (Y-'^) 


6.82 


17.1 


68.2 


'^drag,em (Y-I^) 


812. 


6.99 X 10^ 


2.13 X 10^ 


Tdrag (jr) 


6.76 


17.1 


68.2 


ojt (rad s~^) 


4.70 X 10^ 


4.76 X 10^ 


1.49 X 10^ 


{uJn2/uJTf for to ~ lOVr 


0.209 


0.522 


2.09 


{Urad/uJTf for 6 = 0° 


7.59 X 10^3 


2.65 


1.79 X 10^ 


{uJradluJT? for 6 = 60° 


1.42 X 10-2 


1.31 


1.38 X 10^ 



Note. — We assume the grain parameters of Table ||, p = 3 gem ^, ai = 1.745, 6 = 2, 
Tin = lO^cm-3, T = 45K, / = 1/3, n(H)/nH = .01, E^^ = 0.2eV, / = 10 A, = 45K, 
Urad = 240misrf, and 7 = 0.3. 
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Fig. 1. — Two examples of macroscopic targets, (a) A target on which an isotropic radiation 
field would exert a torque along axis ai. The target is perfectly refiecting except for the 
regions painted black, (b) A perfectly-reflecting target on which an isotropic radiation field 
would exert zero torque, radiation antiparallel to ai would exert a positive torque parallel 
to ai, and radiation parallel to ai would exert a torque antiparallel to ai. 
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Fig. 2. — Grain orientation in "scattering coordinates' 
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Fig. 4. — Extinction and polarization efficiency factors for the grain shape shown in 
Fig. ^ with aeff = 0.2/im, spinning around the principal axis ai, with radiation propagating 
perpendicular to ai. The polarization cross section peaks at A = 0.5/im, with Qpoi/Qext = 
0.086 for this particular grain. 
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Fig. 5. — Same as Fig. |, but for agg = 0.05/im, with Qpoi peaking at A = 0.2/xm, with 

Qpol/ Qext = 0.15. 
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A (/xm) 

Fig. 6. — The component of the radiation pressure efficiency vector Qpr{Q, 0) along k, the 
propagation direction of the incident radiation, as a function of wavelength A, for several 
values of the angle between k and the grain axis ai, and for two different grain sizes, 
ttcs = 0.2/im and 0.05/xm. Qp,. is averaged over the angle f3 measuring rotations of the grain 
around ai. The peak near 10/xm is the silicate absorption feature. 
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Fig. 7. — Three components of AQr for O = 30°, as a function of A, for acs = 0.2/im 
grain. Results are shown for DDA calculations using N = 6656, 22464, and 53248 dipoles. 
It is seen that the N — 22464 dipole array provides a good approximation to the target for 
A ^ 0.15//m. 
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0.1 1 10 

A (/xm) 

Fig. 8. — AQr(©, 0) -ai, where Qr is the radiation torque efficiency vector and ai is the 
principal axis with the largest moment of inertia, for — 0.2/xm grain and various values 
of the angle 6 between ai and k. 



-40- 




-41 - 



0.04 



0.02 



<a3 

A 

V 







0.02 



0.04 



0.06 




agfj^0.2/xm 



UV 

7500K 
4000K 
3000K 
ISRF 



-0.5 





cos© 



0.5 



Fig. 10. — The component of the spectrum-averaged radiation torque efficiency vector (Qr) 
along ai for a^s = 0.2/im grain, multiphed by A, as a function of cos0, for the five radiation 
fields of Table |^. (Qr) is averaged over rotations of the grain around ai. 
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Fig. 11. — Same as Fig. |To]but for = 0.05/im. The results for the "UV" radiation field 
have been multiplied by a factor 0.1. 
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Fig. 12. — Same as Fig.^ except showing the component of the spectrum-averaged radiation 
torque efficiency vector along 63. This torque would tend to cause precession of J around k. 
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Fig. 13. — Same as Fig.|12| except for a^s = 0.05/im 
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Fig. 14. — Same as Fig.|TO| except showing the component of the spectrum-averaged torque 
efficiency vector (Qr) along ai x 63, for a^s = 0.2/im grain. If ai || J, this torque (if positive) 
would tend to cause alignment of J with k. 
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Fig. 15. — Same as Fig.|Tj but for Oefr = 0.05/im. The results for the "UV" radiation field 
have been multiplied by 0.1 . 



